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The Hamiltonian treatment for the collapse of thin shells for a family of Lanczos- Lovelock theories 
is studied. This formalism allows us to carry out a concise analysis of these theories. It is found 
that the black holes solution can be created by collapsing a thin shell. Naked singularities cannot 
be formed by this mechanism. Among the different Lanczos-Lovelock's theories, the Chern-Simons' 
theory corresponds to an exceptional case, because naked singularities can emerge from the collapse 
of a thin shell. This kind of theory does not possess a gravitational self-interaction analogous to the 
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Newtonian case. 
PACS numbers: 



In a recent work Q], black hole solutions in a particular class of Lovelock's gravitation theories were studied. These 
■ theories were selected by requiring that they have a unique Anti-de Sitter (AdS) vacuum with a fixed cosmological 
constant. This strongly restricts the coefficients in the Lanczos-Lovelock (LL) action 0. For a given dimension d, 
the La grangians under consideration are labelled by an integer k = 1, [^rp] 1 , where the Einstein-Hilbert (EH) 
Lagrangian corresponds to the case k = 1. For k = [^f±], we must distinguish between even and odd dimensions, 
I ■ because the theories are different. When d is odd, the corresponding Lagrangian is given by the Euler-Chern-Simons 
form (CS) for the AdS group, whose exterior derivative is proportional to the Euler density in 2n dimensions 
For d even, the Lagrangian reads as the Born-Infeld form (BI). In this case the expression for the (LL) action is 
proportional to the Pfaffian of the 2-form R ab and, in this sense, it has a Born-Infeld- like form These two cases 
are exceptional, because they are the only ones which allow sectors with non-trivial torsion . For d > 7 there exist 
other interesting possibilities, which are different from EH, BI and CS. For example, the theory with k — 2 has been 



D \ studied by several authors in different scenarios [3, la 111 III 

. In this LL theories for any dimensions and k, there exist well-behaved black hole solutions. However, we must 
differentiate between cases with odd and even k, because in theories with even k an, additional solution appears, 
• i-h , which represents a naked singularity. 

It is interesting to study the black holes formations through gravitational collapses of thin shells. In the usual thin 
shell treatment [H], IeH , the analysis of collapse is based on the discontinuities of the extrinsic curvature of the 
world tube of the collapsing matter. However, the implementation of the Israel formalism in this kind of theories (LL) 
is very difficult, because the action contains high powers in the curvature and, therefore, in the extrinsic curvature. 
In this formalism the complicated analysis of collapse makes the treatment quite unattractive. 

Another approach in studying matter collapses is the Oppenhcimer-Snyder formalism, which was applied by Ilha 
ct al [Til Hl| to the case of a homogeneous collapsing dust, where the inner metric is described by the Friedmann- 
Robertson- Walker line element, and the external metric corresponds to the solution of the fields equations in BI or 
CS theories. In Ref. 0], the authors discussed the formation of a naked singulary in the CS theory. 

On the other hand, an alternative way to study gravitational collapse of thin shell is the Hamiltonian treatment. 
This treatment was applied in Ref. 0] to the Einstein-Hilbert gravity, where the direct integration of the canonical 
constraints reproduces the standard shell dynamics for a number of known cases. In particular, it was applied in detail 
to three dimensional spacetime and the properties of the (2+l)-dimensional charged black hole collapse was further 
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1 Here [x] is the integer of x. 
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elucidated. The Hamiltonian treatement was also extended to deal with rotating solutions in three dimensions. The 
general form of the equations of the shell dynamics implies the stability of black holes. As far as, black hole in this 
model cannot be converted into naked singularities by any shell collapse processes. 

In this work we will extend the Hamiltonian formalism to our approach to the theory of LL in high dimensions, 
and particularly to the theories described in Ref. [lj . This formalism permits us to analyze the black holes formations 
in an economical way. 

The plan of the paper is as follows. Section II briefly reviews the LL action and its spherically symmetric solution. 
Section III analyzes the collapse of a spherically symmetric shell under the Hamiltonian formalism. Finally, section 
IV is devoted to conclusions. 



II. LL ACTION 



The LL action is a polynomial of degree [d/2] in the curvature, which can be also written in terms of the Riemann 
curvature R ab = duj ab + ui®uj cb and the vielbein e a as 2 



I (g) = K 



[d/2] 
p=0 



.R 



(1) 



where a p are arbitrary constants. In the first-order formalism, the action (|TJ is regarded as a functional of the vielbein 
and the spin connection, and the corresponding field equations obtained varying with respect to e a and ui ab reads 

[^] 

(d-2p)a p e abl ... bd _ 1 R b ^...R b ^-^e b2 "+\..e b ^ = 0, (2) 

p=0 



r d- 1 1 



E 

p=0 



p{d- 2p) a p e 



p£abc 3 ...c d 



R C3C4 ...R C2 " 



2p r pC 2p +l e C2p+2 



= 0. 



(3) 



Here T a = de a + u; b e b is the torsion 2-form. 

Note that in even dimensions, the term L^ d l 2S > is the Euler density and, therefore, does not contribute to the field 
equations. However, the presence of this term in the action -with a fixed weight factor- guarantees the existence of 
a well-defined variational principle for asymptotically locally AdS spacetimes [jjl U_\ ■ 

The first two terms in the LL action JJJ are the cosmological and kinetic terms of the EH action. Therefore, 
General Relativity is contained in the LL theory as a particular case. 

The linearized approximation of the LL and EH actions around a flat, torsionless background are classically 
equivalent |19|. However, beyond the perturbation theory the presence of higher powers of curvature in the Lagrangian 
make both theories radically different. In particular, black holes and big-bang solutions of J2J have different asymptotic 
behaviors from their EH counterparts. Hence, a generic solution of the LL action cannot be approximated by a 
solution of Einstein's theory. 



A. Static and Spherically Symmetric Solutions 



Consider static and spherically symmetric solutions of equations |J2J and 10. In Schwarzschild-like coordinates, the 
metric can be written as 



ds 2 



-N 2 (r)f(r)dt 2 



dr 2 

_____ I r 2 /7Q 2 



(4) 



2 Wedge product between forms is understood throughout. 
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Replacing this Ansatz in the field equations J5J) and ® leads to the following equations for N(r) and / 2 (r 



dN 
dr 



= 0, (5) 



dr 

Integrating equations JSJ and © yields 



( r^J2(d-2 P )a p ( 1 —f) P ) = 0. (6) 



p 



(7) 

M + C 

K (d-2)!r> d _ 2 

where the constant of integration relates the coordinate time to the proper time of an observer at spatial infinity. 
We will assume it equal to one. The constant M stands for the mass up to an additive constant Co, which is nonzero 
only in the case of CS theories |l|. 

Equations JHJ) corresponds to the solution of field equations, which is a polynomial in f 2 (r), so many roots for / 2 (r) 
with the same mass will exist, but, with different asymptotical behavior. This means that (J2J possesses, in general, 
several solutions with a constant curvature in the asymptotical region, making the value of the cosmological constant 
ambiguous. In fact, the cosmological constant could change in different regions of the same spatial section, or it could 
jump arbitrarily as the system evolves in time |20| . 

These problems are overcome by demanding that the theory have a unique cosmological constant Q. In order to 
satisfy this condition, we choose the coefficients a p 's as follows: 



k 



p < k 



and 



2(d~2y.n d - 2 G k 
With this choice, f 2 (r) adopts the following form: 

r_ 2 , , (2G k M + 5 d _ 2k ^ 1/k 
P 



^{pj'"- , (9) 
,p>k 



(10) 



f(r) = 7a + 1 ~ X ( ILL i ' ) , (ID 



where \ = (±l) fc+1 . For even k, the ambiguity of sign expressed through \ in l|ll|) implies that there are two possible 
solutions, provided M > 0. The solution with \ = 1 describes a black hole with an events horizon surrounding the 
singularity at the origin. The solution with \ = ~ 1 h as a naked singularity with positive mass. If k is odd, there is 
no ambiguity of sign because % = 1. Therefore this solution corresponds to a black hole with positive mass. 

From cq. (|ll|l . it is observed that for k — 1, the AdS black hole solution for EH in <i-dimcnsional is recovered. 
The black hole solutions corresponding to BI and CS theories are obtained also from expression Qllfl. setting 



3 In the first order formalisrn, the field equations imply that the torsion vanishes, except for BI and CS theories, so that, it is not 
necessarily to set T a = (j. However, for static and spherically symmetric configurations the equation implies that the torsion 
must vanish in these cases as well. 
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III. COLLAPSE OF THIN SHELLS 

Let Eg be a time- like hypersurface, which represents the evolution of a thin shell [TllIT^ ]. This hypersurface divides 
the spacetime into two regions; the interior denoted by y' - ) and the exterior denoted by V^ + \ respectively. Each 
of these regions contains Eg as a part of its boundary. We introduce into Eg a set of intrinsic coordinates p a , where 
the Latin indices go from to d — 2, and in the regions V^ - ' and V^ + \ the independent coordinates x°_ and x? are 
introduced, so that the parametric equations for Eg in these charts are x" (p a ) and x" (p a ), respectively. 

At each point on Eg there exists a unit space-like vector normal to Eg and pointing from TA _ ) to V^ + \ and 
d — 1 vectors e" = dx a /dp a tangential to Eg in the directions of the coordinates p a . 

The time-like hypersurface Eg represents the evolution of an infinitesimal d— 2-dimensional matter thin shell. There 
is no matter outside the shell. Therefore, the matter moves only on the shell, so that its d-velocity u a is normal to £ A 
and vanishes outside Eg. Moreover, an observer on the shell can refer the movement of matter to the reference points 
(p 1 , .... p d ~ 2 ) and the reference time p° = r, and thus, the velocity is described by an intrinsic vector u a . The vectors 
u a and u a are joined by the relation u a = e"u a . 

The mechanical properties of matters are described by the surface energy-momentum tensor T^ v , which is normal 
to £ A and it vanishes outside Eg. For an observer on of Eg, the tensor is described by the intrinsic coordinates. 
For an ideal fluid, the intrinsic energy-momentum tensor has this form: 

Tab = <JU a U b - f(h a b + U a U b ), (12) 

where a means the rest surface mass density of the shell and f the surface tension 4 . Since, the tensor T ab is confined 
into the hypersurface Eg, it satisfies the continuity equation T^ b = 0. Multiplying this tensor by u a , we obtain the 
following relation 

K)/« - ru% = 0. (13) 

This equation can be seen as the equation of state of the matter on the hypersurface Eg. 

The next step is to introduce a timelike ADM foliation, E t , of the spacetime. The foliation intersects the world 
tube of the collapsing matter, which corresponds to the thin shell of the E t at the time t. As usual, the metric tensor 
is decomposed by using the basis N» = {N^,N l ), where N x represents the lapse and N l the shift function. In this 
basis the line element in the coordinates x a of the regions y(~) and takes the form 

ds 2 = -{N^fdt 2 + g lj (N l dt + dx l ){N^dt + dx j ). (14) 
In the presence of matter, and since are Lagrange Multipliers, the total Hamiltonian becomes 



where H± and Hi arc defined by 



and 



H = N^Hi. + N*7U, (15) 



n^ = n { l ) +n { r\ (16) 



Hi = Hf +H ( i m \ (17) 

respectively. Here Ti^f 1 and Ti.^ correspond to the Hamiltonian terms related to the gravitational field of the LL 
action [2l]], and are given by 

n-l(s) _ „ \ " (j ~ 2p)!<3p x i 1 ...i 2p p iij 2 p hp-lhp ho\ 



4 Here f denotes the surface tension and r, the proper time. 
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and 

^ = -274., (19) 

where tt^ are the conjugate momenta to the metric tensor of the intrinsic tensor metric g, its determinant; and 
Rijki corresponds to the tensor curvature, which can be written in terms of the geometric quantities of E t as 

Rijki = Rijki + K lk Kji - KuKjk, (20) 

with Rijki stands for the curvature tensor of the leaf E t of the foliation and Kij is the extrinsic curvature. 
The momenta are defined in terms of extrinsic curvature — j^i^i/j + ^j/i ~~ 9ij) as 

tt* - r. /77 V ^ d - 2p ^- a P \ " n , f»i-»2«-i2„-i p hh r> ha-ih, K h>+i K h v -i f 2 i\ 

p s— 

where 

(_4)P" S 



D s( P ) = 



«![2(p -«)-!]!!' 



The matter components and H -"^ , are given by 



= (22) 



U^ l) = 2VgT ±i , (23) 

where _L corresponds to a contraction with a normal vector to the hypersurface E t , = (— N- 1 , 0, 0). 

In what follows we restricts ourselves to the spherical coordinates. We will use the proper time r and spherical 
angles as intrinsic coordinates of the hypersurface E^; p a = (r, 6 1 , 9 d ~ 2 ). The motion of the shell is expressed 
by the equation r = R(t). The derivative with respect to r is denoted by a dot. The line element of in this 
coordinates is expressed by 

dsl = -dr 2 + R 2 { T )dtt 2 d _ 2 . (24) 
The interior and exterior line element with spherical symmetry are given by 

ds 2 _ = -f_{r)dt 2 _ + fl 2 (r)dr 2 + r 2 dtt 2 d _ 2 , r < R(t), (25) 

and 

ds% = -fl(r)dt 2 + + f+ 2 (r)dr 2 + r 2 dQ 2 _ 2 , r > R(r). (26) 

Interior and exterior coordinates match continuously on the Ej, but it is found that t- ^ t+. In these coordinates 
the vectors u a and £ a are given by 

u a = (-^,i?,0,...,0), (27) 

r - (p7,0,.,0), (28) 
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where 

7 = \jp + R 2 . (29) 
Substituting (|25|l and (|26() into the Hamiltonian generator H± l|16[) . we obtain 

^ = - K ~^^i X> - 2 ^K ( i 7 / ! ) P } + 00) 

We are interested in integrating out the constraint H± = across a radial infinitesimal length centered in the shell 
position r = R(t) to a constant time. In this form, it is possible to express the discontinuities of geometry in terms 
of the projected stress Tj_x< It is easy to prove that Tj_j_ has the same form that the one obtained in 0], due to the 
symmetry of the thin shell. Finally T±± is given by 

Tj_ ± = <ry<5 (r - R(t)) . (31) 
Integrating in the radial direction, from R + e and R — e, we obtain in the limit e — > 

«(d- 2)!*5> - 2p)a p | (^J^)" - ( 1 ^ M ) P } - ^(7 + + 7-)- (32) 

This expression has been seen as the generalization of the equation obtained for GR. 
From expression J^J, it is found that 

M+-M_ = im( 7+ + 7-), (33) 

where m = Qd-2R d 2 <^- Expression 1|33[) is the same to that obtained from the GR case If a > then M + > M_, 
this means that if M_ is the mass of a black hole inside of the shell, the final mass of the black hole will be greater, 
therefore its events horizon increases. 

In order to complete the present picture of a radial collapse, it is necessary to analyze the consistency of the 
remaining nonvanishing components of the Hamiltonian treatment related to the radial and angular components. The 
angular contribution of the constraint (|17|) are identical to zero. Because the radial component is not identically zero, 
it is necessary to evaluate 

H r = -2rf r/j + 2^/gT ±r , (34) 

which yields 

U r = ^(d-2y.± ^~_ 2 ^ ^rD^f^ (1 - f)* ± {r ^-\arf^ )-^^aS(r-R(r)), 

p=0 ^ " s=0 •* 

where T corresponds to the determinant of the angular metric. One can expect that H. r to be proportional to %x, 
since l|16|) already provides the equation of motion for R(t) . It would be interesting to explicitly see that this indeed 
occurs. But due to equation (|35J) . in the general case it is complicated (it is not possible to integrate) to perform this 
point. However, it is straightforward to prove that the correct Einstein-Hilbert limit is obtained when k = 1 |l6j |. 

The acceleration of the thin shell is obtained from equation ll-i-it by differentiating with respect to proper time r, 
which yields 

= -^TT7^T-, ( - (d - 2) ^_ 2 i? d - 3 f 7+7 _, (36) 



2{M+-M_) V dR '~ dR 
Notice that we need the explicit forms of /_ and /?. The form of /_ will be 
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f-(r) = ¥ + l. (37) 

The form of /? must be split into two cases; the case d — 2k — 1 ^ and the case d — 2k — 1 = 0. In the latter case 
there exists a gap in the mass, in which the vacuum corresponds to M_ = — (2Gfc) _1 . 

• If d — 2k — 1 7^ 0, with M_ = and M + = M, then equation l|33|l takes the form M = hm (7+ +7-), so that if 
a > 0, then M > 0. The acceleration is given by 



mR = --p?R- (d- 2)n d _ 2 R 
V 

X (d - 2k - 1) m 2 



d-3 



2k 



T 7+7 _ 
2Gk 



M k ' 1 R d - k - 1 



i/fe 



7- 



(38) 



The first two terms of l|38|l correspond to the acceleration due to AdS, and the interaction of tangent tension 
on the thin shell. If \ = 1, for a black hole solution, then R < 0. Therefore, in this way the thin shell always 
collapses to a black hole. 

On the other hand, when x = — 1 we might think that a naked singularity , however due to that the latter 
term a Eq. I|38|) is positive and thus a repulsive gravitational force appears. It could be shown that this force 
therefore domines over the other terms when R — > 0. Therefore, naked singularity cannot be formed through a 
thin shell collapse. In order to see this point we consider 



M = -( 7+ + 7 _), 

where 

7± = \Jr 2 + p±, 

and m = fld-2R d 2(7 with M > 0. In this case and f+ are given by 



f 2 = 1 + — 



rf 

p 



and 



R 2 ( 2G k M 



p \ R d-2k-l 



l/k 



respectively. For a naked singularity it is necessary that /? > f 2 : > 0. 



From (1231 we obtain i? 2 



(39) 



(40) 



(41) 



(42) 



R 2 



M (2 1 - 2k G k M 1 - 



k \ i/fe' 



/ 2 



1 



(43) 



from which we could read an effective potential. 

In order to give either a quantitative and qualitative discussion of this potential, let us write Eq. I|43|l in the 
form 



where 



R = a + a- 



_ M (2 1 ~ 2k G k M 1 



k\ Vfe 



,R 2 

±\l-w + i. 



(44) 



(45) 
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R 



FIG. 1: This plot shows behavior of our effective potential and their turning point. 



The physical regions are defined by R 2 > 0, which implies that both a + and a_ have the same sign. Note that 
for the case under study we have 7± > 1, due to f± > 1, this means that the sum 7+ + 7- > 2. Therefore, 
from Eq. i|39|) we obtain M > m. In order to characterize the physical regions, we need the behaviors of the a± 
parameters. We will simplify our study to the dust case which means m = Constant. Our results are shown in 
Fig.l from which we could see different regions, 



I. 


R 2 


> 


0, for R < Ri, 


II. 


R 2 


< 


0, for R 1 <R<R 2 


III. 


R 2 


> 


0, for R 2 < R< i? 3 


IV. 


R 2 


< 


0, for R 3 < R, 



Note that the classical allowed regions are I and III. In region I note that when R — > R 2 — * 00. In this case, a 
naked singularity emerge from a collapse of a thin shell. However we could prove that the motion in this region 
is prohibited by causality. In fact, from Eq. I|52fl ( see appendix), we obtain 



M {2 1 ~ 2k G k M 1 



\-lkn n/ri-k \ 1 / k 



7- = --"*( Rd -2*.-r ) (46) 



and the condition 7_ > 0, implies 



m >m V R ^l x I'''- (47) 



M / , 2 1 - 2k G h M 1 ~ k 



This mean that it must exist a minimum radius, which we denote by R = i?* in order not to violate the 
conservation law expressed by Eq. (|39|l . It is direct to prove that i?i < i?* < R2, which implies that a + > 0. 
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Therefore, the thin shell "movement" is allowed only in region III, this is R2 < R < R3. So, naked singularity, 
in the case of dust, cannot be formed through a thin shell collapse. 

For the general case, is not simple to study equation (|43|l . since it requires to solve equation l|57|l (see appendix). 
Also we should know the relation between a and f (the state equation for the shell), in order to obtain m = m(R). 
However, we may argue that equation (|4T?)l is of general character, so also it applies for the general case. The 
a minimum radius must exist, so that the thin shell does not violate equation <|4(>p . Therefore, it must exist a 
turning point. So, naked singularity, in the general case, cannot be formed through a thin shell collapse. 

It is easy to see that in the limit R — *■ the first term vanishes. If we consider the particular case of dust, i.e. 
f = 0, eq. H57fl implies that m = Const., therefore the last term goes to infinity for R — ► 0. Thus, for the dust 
case the thin shell does not collapse to R = 0, because the acceleration becomes very strong. 

At this point, if we considered Einstein- Hilbcrt limit (k — 1) in equation (|38[l . the last term would be reduced to 
-(d— 3)Gm 2 /R d ~ 2 , thus corresponding to the Newton gravitational interaction. Therefore, for k ^ 1 this term 
will be a generalization of the Newton gravitational interaction, with an effective gravitational constant given 
by 



1 / 2G k x 1/fe 



k \M k ~ 1 



(48) 



• If d — 2k — 1 = 0, it corresponds to the CS theory with M_ = — (2Gfc)~ 1 . In this case, expression IjHHJ) takes 
the form M + (2G k y 1 = \m (7+ + 7_). If a > 0, then it is implied that M > -(2G fe ) _1 ; therefore, the naked 
singularities with negative mass can emerge from the collapse of a thin shell. Here, acceleration is given by 

mR = {d-2)n d -2R d ~ 3 Tj +1 -, (49) 

where R < 0, and thus the thin shell always collapses. It is observed from eq. (|49|l that a term analogue to the 
Newton gravitational interaction does not appear. 

IV. CONCLUSION AND REMARKS 

We have developed the Hamiltonian formalism for the collapse of thin shells in Lanczos-Lovelock theories, and we 
presented given a concise analysis of the theories described in Ref. [lj . We show in these theories that the black holes 
solution can be created by collapsing a thin shell and naked singularities cannot be formed by this mechanism. On the 
other hand, if we consider theories with k ^ 1, these exhibit a generalization of the Newton gravitational interaction, 
and effective gravitational constant becomes given by 



1 / 2G k \ 1/k 



k \M k - 



(50) 



Also we have shown that when we take the Einstein-Hilbert limit (k = 1) in equation (|38[1 . the last term is reduced 
to -(d — 3)Gm 2 / R d ~ 2 , which coincides with the Newton gravitational interaction. 

Nevertheless among the different Lanczos-Lovelock's theories, the Chern-Simons theory exhibits an exceptional 
behavior, since naked singularities can emerge from the collapsing of a thin shell. This kind of theory does not possess 
a gravitational self-interaction analogous to the Newtonian case. 

It is straightforward to prove that in the case of electrically charged thin shells, the general form of eq. %\2\ . 
governing the radial collapse in d dimensions, remains the same because the electromagnetic stress tensor contributes 
with a finite jump value on the Tl± and, therefore, does not contribute to the radial integral of H±. Moreover, when 
we consider the charged case in CS theories, a mechanism that prevents the naked singularities formation appears. 

Finally, we conjecture, by virtue of the results from Ref. |16| that the presence of an angular moment in 2 + 1 
dimensions prevents the formation of naked singularities; thus, in higher dimensions, the angular moment could 
prevent naked singularities formations in CS theories. 
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VI. APPENDIX 



From Eq. (15). i.e. 



R 2 



M 



m 



2 1 ~ 2k G k M 1 



R 2 
P 



(51) 



we will obtain the acceleration of the thin shell given by Eq. (|38p. 
Let star from 



M 

7± = — ± 



and 



2 1 - 2k G k M 1 - k ^ 1/k 

R d-2k-l 



)/« - ru% = 0. 



(52) 



(53) 



This equation can be seen as the equation of state of the matter on the hypersurface . Using spherical coordinates 
and the identity T? = d a (^/—hT a )/^/—h, where h is the proper metric determinant, we obtain 



The proper metric is given by 



ds 2 



-dr 2 + R 2 {T)dn 2 d _ 2 



(54) 



(55) 



with y/—h = R d 2 V / T, where T is the angular metric determinant. Besides we consider radial collapse then u a 
(1, 0, 0), moreover consider proper coordinates given by p a = (t, 9 1 , 9 d ~ 2 ), then we get 



dr 



R d - 2 Vfa) =(d- 2)R d ~ 3 r Vf R, 



where R — dR / dr and since dVT/dr = we have 

d 



dr 



{R d ~ 2 a) ={d- 2)R d - 3 TR 



(56) 



(57) 



Multiplying this latter equation by the angular volume of the unit sphere, fld-2 (for d — A Q 2 — 47r) and defining 
m = fld-2R d ~ 2 cr, we obtain 



^ = (d-2)R d - 3 rR. 
dr 

Due to, symmetry we have m = m(R) and d/dt — Rd/dR, then 

dm 



dR 



{d-2)R d - i T 



(58) 



(59) 
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From which we get 
2RR = 



2 • 
p RR 



M 



2 1 - 2k G k M 1 - k ^ 1/k ' 

ud-2k-l 



M dm { 2 1 - 2k G k M 1 - k \ 1/k dm t _(d-2fc-l) (2 1 - 2k G k M 1 - k N 1/fc 



to 2 dr V R d - 2k - x J dr 
For the (— ) sign, we obtain the following 

2RR = — —RR 

l z 

M dm f2 1 ~ 2k G k M 1 - k ^ 1/k 



R 



+2-t- 



m 2 dr 



ftd-2k-l 



dm 
~dfr 



m- 



(d-2*-l) (2 1 - 2k G k M 1 ~ k ^ 1/k 
R d-k-\ 



R 



we could rewritten this equation in the form 
1 



RR 



-RR 



_1_ 
m 



M 
m 



l 1 ~ 2k G k M 1 ~ k x 1/k 



Jjd-2k-l 



dm (d-2k- 1) (2 1 - %k G k M 1 -' 



dr 



R d - 



k-l 



1/k 



mi?7_. 



For (+) sign, we obtain 



RR - 1 RR 7 - 7 + dm 
JXJX — —-77 lire — 



m dr 

■\\-ikn A/fl-fe \ x l k 



(d-2k-l) (2 1 ~ 2k G k M 1 



d-k-l 



but dm/dr — Rdm/dR and using Eq. (|5 

I? 

multiplying by m, we obtain finally 



~R- {d-2)R d -H^2± 
V m 

(d-2k- 1) {2 1 - 2k G k M 1 ~ k 

Rd-^ 1 



i/fc 



?717_ 



mi? = ~R- (d-2)R d - 3 Tj-^ 



(d-2k- 1) (2 x - 2k G k M x - k 



■fc-i 



m 7_ 



which corresponds to Eq. Ij38(l in the main text. 
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